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O |' Abstract 

<**> I We present the analytical calculation in perturbative Quantum Chromodynamics of 

the a^ contribution to the Ellis-Jaffe sum rule for the structure function g\ of polarized 
deep inelastic lepton-nucleon scattering. 



PACS numbers: 12.38.Bx, 12.38.-t, 13.60.Hb, 13.85.Hd 



Measurements of the polarized nucleon structure function g\ during the last 20 years have 
revealed an internal spin structure of the nucleon that is surprisingly different from constituent 
quark model expectations. The discovery of the disagreement between the combined EMC- 
SLAC data (l|, §] and the constituent quark model expectation || J dxgi(x,Q 2 ) « 0. 1 5 1 (7^4 1 
attracted a lot of attention and triggered intensive research in the field of polarized deep 
inelastic scattering. More recently, deuterium scattering data of the SMC Q and of the E143 
H collaborations, and 3 He scattering data of the E142 || and E154 Q collaborations also 
allowed the determination of the neutron sum rule J dxgi(x, Q 2 ) and the Bjorken sum rule 
H /o dx[g\(x, Q 2 ) — gf(x, Q 2 )]- At present the Bjorken sum rule is confirmed at the 8% level 
P] which is an important experimental test of Quantum Chromodynamics. 

Higher order perturbative QCD corrections to the sum rules are crucial for an accurate 
and reliable confrontation of these sum rules with experimental data, see e.g. Ref. [10|. For 



the Bjorken sum rule, the a s correction [11], the a 2 correction |12|| , and the a 3 correction [13] 
have been calculated in the leading twist approximation. Higher twist corrections have also 
been calculated [14|. The Ellis-Jaffe sum rule J Q dxg\{x,Q 2 ) for the proton and neutron 



was calculated to order a s [15] and to order a s JT6| in the leading twist approximation. Power 
corrections were calculated in [jl7| . 

In this article we obtain the order a^ contribution to the Ellis-Jaffe sum rule in the leading 
twist approximation for massless quarks. We perform the calculations using dimensional 



regularization [18] in D = 4 — 2e space-time dimensions and use the standard modification of 
the minimal subtraction scheme fl^| , the MS-scheme ] p0[ . 

Polarized deep inelastic electron-nucleon scattering is described by the hadronic tensor 

WV = ^ J ^ze^ z (p,s\J^z)J v {Q)\p,s) 

= [~9iw + -^J Fi(x,Q 2 ) + (pim r qA \Pu 2~*0 — -^(^Q 2 ) 

( Sct 1 r&\ , s a p ■ q - p a q ■ s 2 \ 

+ie»v P eq P l^—gi(x, Q ) + ^—^ g 2 (x, Q )j (1) 

Here J M = J2i£i Zi^il^i ls the electromagnetic quark current where ej = 2/3, —1/3, —1/3, • • • 
is the electromagnetic charge of a quark with the corresponding flavour, x = Q 2 /{2p- q) is the 
Bjorken scaling variable and Q 2 = —q 2 is the square of the transferred momentum. \p, s) is 
the nucleon state that is normalized as (p, s\p', s') = 2p° (2ir) 3 5^' (p — p')5 ss i. The polarization 
vector of the nucleon is expressed as s a = U(p, s)^ a ^U{p, s) where U(p, s) is the nucleon 
spinor U(p, s)U(p, s) = 2M. 

In the present article we will focus on the first Mellin moment of the structure function 
g\. Moments of deep inelastic structure functions can be expressed [^IJ (for reviews see Refs. 
[22, ^3|) in terms of quantities that appear in the operator product expansion (OPE) of the 
two currents J^. For the first moment of the structure function g\ we need to consider the 
following expression for the OPE of two electromagnetic currents 



, I dze^ z T{J^z)J v {G)} Q2 =°°e^ p(T |§ 



+ C s (log(^),a s (v 2 ))4(0) 



+ ••• (2) 



where everything is assumed to be renormalized (with \x being the renormalization scale). 
Jp ,a (x) = il)^ a ^t a "ijj{x) is the non-singlet axial current, where t a is a generator of the flavour 
group, and J„(x) = J2i=i V'jTo-TsV'j^) is the singlet axial current. We use the notation 

as{l " ] = -m^ = -*r (3) 

for the QCD strong coupling constant, where g{^ 2 ) is the coupling constant of the QCD 
Lagrangian. 

The dots on the r.h.s of Eq. (||) indicate higher spin and higher twist contributions that 
do not contribute to the first moment of g\ in the leading twist approximation. The strict 
method of the OPE ensures [24] that the OPE of two gauge- invariant currents can only 



contain gauge invariant operators with their renormalization basis. For example, the known 
twist-two and spin-one axial gluon current K a = Ae apLV p{A a d u A a + gf abc A a A b v A i ]/3) also has 
the necessary quantum numbers, but it cannot contribute to the above OPE because it is 
not gauge invariant. Since the formalism of the OPE gives the strict theoretical basis to the 
parton model with the parton distributions being the matrix elements of the corresponding 
operators, there is strictly no well founded way to get the polarized gluon distribution into 
the considered sum rule. 

Since the non-singlet coefficient function C a depends trivially on the flavour number a (see 
e.g. Ref. |25|) it is possible to factorize a flavour number independent coefficient function for 
which we use the standard notation C ns . 

The Ellis-Jaffe sum-rule is expressed as 



f'dxg^^Q 2 ) S™" - <^~ ' ' 

Jo 



C ns (l,a s (^))(±-|< M | + -a 8 ) 
12 6b 



+C s (WQ 2 ))exp / da',l±jd. -a (^ 2 ) (4) 

\Ja a (u2) B(a') J 9 



3 (M 2 ) 

where the plus (minus) sign before \gA\ corresponds to the proton (neutron) target. The 
proton matrix elements of the axial currents are defined as 

\9A\sa = 2{p,s\Ja' 3 \p,s) = (Au- Ad)scr, 

a 8 s a = 2y/3(p,s\J^ 8 \p,s) = (Am + Ad- 2As)s a , (5) 

ao(v?)s a = (p, s\J%\p, s) = (Au + Ad + As)s a = AT,(fi 2 )s a . 

Here \gA\ is the absolute value of the constant of the neutron beta-decay, gAJgv = —1.2601 ± 
0.0025 1 26]. a% = 0.579 ±0.025 [^, ^7]] is the constant of hyperon decays. The matrix element 



of the singlet axial current ao(/U 2 ) will be redefined in a proper invariant way as a constant 
So below in Eq. (|l0|) . 

We use the notation /S.q([i 2 )s a = {p,s\q'y a 'y5q\p,s), q = u,d,s, for the polarized quark 
distributions. We omit the contributions of the nucleon matrix elements for quarks heavier 
than the s-quark but it is straightforward to include them. We also avoid to introduce the 
contribution of the polarized gluon distribution Ag in the expression for the matrix element 
of the singlet axial current ao(/x 2 ) for the reason given below Eq. (|3]). 

P(a s ) is the beta function that determines the renormalization scale dependence of the 



renormalized coupling constant. It is presently known at four loops [28] in the MS scheme 

(JO 

P(a s ) = j^ = -Poa 2 s -(3 1 a 3 s -p 2 a 4 s -l3 3 al + 0(a 6 s ) 

with the SU(3) values 

fa = 11 - §n/ 
/3i = 102 - fn f 

a _ 2857 5033 ^ , 325^2 

a _ ( 149753 i qcc/iA A ( 1078361 , 6508 a 
P3 - { — h 3564C 3 J - ^ leg + ~TT^ 



»/ 



i / 50065 i 6472 a \ „2 , 1093 ^3 lp\ 

+ \~v&T + -8I-C3J n f + - w n / (b) 

in which £ is the Riemann zeta-function (£3 = 1.202056903 • • •). 

The singlet anomalous dimension 7 s that appears in Eq. (Q) determines the renormaliza- 
tion scale dependence of the axial singlet current J^, 

'' 'Jl]R = lVlU (7) 



din// 2 

where subscript R means that a current is renormalized. 

The axial singlet current is not conserved - the axial anomaly p9| - and this causes the 
singlet anomalous dimension to be non-zero starting from the order a 2 . 

j s (a s )= 7l a 2 + l2 a 3 s +^ 3 ai + 0(a 5 s ). (8) 

That is why ao(// 2 ) in Eq. (Q) depends on the renormalization point fi 2 and it is therefore 
not a physical quantity. 

The non-singlet axial current is conserved in the limit of massless quarks and the anoma- 
lous dimension for the non-singlet axial current therefore vanishes. This is why the non-singlet 
contribution to the Ellis-Jaffe sum rule Eq. (Q) does not involve an exponential factor and 
that is why the non-singlet matrix elements gA and ag are renormalization group invariant 
(/x 2 independent). 

The non-singlet contribution to the Ellis-Jaffe sum rule is known in the order a? s from [ |l3fl 
where the polarized Bjorken sum rule L dx(gi — g™) was calculated in this order. In this 
article we obtain the order a? s contribution to the singlet part of the Ellis Jaffe sum rule: 



C s (l,a s (Q 2 ))exp 



C s (l,a s (Q 2 )) 



da s —. ,, 

a s (M 2 ) P( a s) 



g«0(/^ 2 ) 



l + a s (Q 



2>>~7l 

/9b 



2\\2' 



+ MQ 2 )) 



-72/3q + liPi + 7i 
2/3 2 



2 ^3 -2 73 /3 2 + 2 7l /? 2 /3 + 2/? l72 /3 - 2 7l /3 2 + 3 7 i 72 /3 - 3 7l 2 /?i - 7l 3 



+(a s (Q 2 )) 
where we introduce the notation 

ao = exp 
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aoO 2 ) 



a (9) 



(10) 



for the renormalization group invariant (i.e. /x 2 independent) nucleon matrix element of the 
singlet axial current. Since ao is the renormalization group invariant it should be considered 
as a physical constant on the same ground as the constants qa and as. That is why from now 
on we will only use the notation ao for the matrix element of the singlet axial current in the 
expression for the sum rule. 

To obtain this singlet contribution, we need to calculate the 3- loop contribution to the 
singlet coefficient function C s and the 4-loop contribution to the anomalous dimension 7 S . 
For the treatment of the 7 5 matrix in dimensional regularization we use a technique described 



in [30 1 which is based on the original definition of 7 5 in [18]. 

In the MS scheme the proper normalization of the axial singlet current requires the in- 
troduction of a finite renormalization constant Z| in addition to the standard ultraviolet 
renormalization constantj^j Z^ s that contains only poles in the regularization parameter e 



Hi 



R 



<7S ys 



[jI\b 






-fiupcr 



^Blvlpla^B 



(11) 



where B denotes the bare, unrenormalized quantity. 

This extra renormalization constant is introduced to keep the exact 1-loop Adler-Bardeen 
form |3l]] for the operator anomaly equation d^J^ = a s ^-z i i V \pG a iliV G a x within dimensional 
regularization. Here G® v = d^A® — d v A" + igf abc A b ^Al is the QCD field strength tensor. 

This finite renormalization constant Z\ affects the results for the coefficient function and 
anomalous dimension but the sum rule, i.e. the combination in the r.h.s. of Eq. (0) is 
independent of the choice of Z\. This fact is evident since the sum rule as a physical object 
can not depend of the choice of the normalization of a non-physical object such as a flavour 
singlet axial current. More precisely, Z\ enters in the coefficient function and anomalous 
dimension as 



7 s 



C s /Z; 



5- 



7 s + f3(a s )—log(ZI) 



(12) 



Please note that ultraviolet renormalization constants coincide in the MS and MS schemes and we can 
therefore use the notation Zms instead of Zj^. 



where C s and 7 s are calculated in the MS scheme without a factor Z\ (i.e. with Z\ = 1). 
Please notice that Z\ enters the anomalous dimension only at the a 2 order, since P(a s ) starts 
from a 2 . Presently the singlet constant Z\ is unknown in the order a^. We will therefore 
obtain C s and 7 s which is sufficient to obtain the sum rule. 

To obtain the 3- loop coefficient function C s we apply the method of projectors [32] which 
gives us the following formula 

^ = ^^%s* /^e^(0|T{^) 7[At7i ,7 P ]9pV'(-p)^^)«/,(0)}|0)|;X tatCd (13) 

Z4 Z MS J 

where 7[ M 7^7p]<7p = lyHvi - S^fi + Qfilu - ft/7/x- The R^g operation performs the renormal- 
ization of the ultraviolet divergences. The infrared divergences that are produced by putting 
p = are removed by the ultraviolet renormalization factor of the axial singlet current Zy^g. 
In this way we need to evaluate 3-loop massless propagator diagrams which is done with the 
package MINCER written for the symbolic manipulation program FORM [J33J 1. The set of 
3-loop flavour singlet diagrams that contribute to this coefficient function is identical to the 



set "q7q7" that previously appeared in Ref. [34]. In this way we obtained the following result 
for C^ 



C s = l + a s C F {-7) 
,3 



2 'C A C F (-^) + ^n / TH^ + 16C 3 ) + C|(|) 



+a„ 



^ ^ ^,41215 1124,., _ ^ 2 ,13004 656 _ 
C A C F n f T F (-^- + — Cs) + C A C 2 F {— — C 3 ) 

., , 13021 > 440 A , „ 2rTl2/ 9404 256 

.-• , r f m aonM , nSf 1397 



+C\C F { — + 56C 3 + — Cs) + C F n}T F {-— — Cs) 



+C F n f T F (-— - 320C5) + C F { — + 96C 3 )j (14) 

where C F = 4/3 and Ca = 3 are the quadratic Casimir operators of the fundamental and 
adjoint representation of the colour group SU(3), T F = 1/2 is the trace normalization of the 
fundamental representation and n/ is the number of (active) quark flavours. The Riemann 
zeta function is written as C, n . It is interesting to mention that diagrams with one external 
photon in a closed quark loop have colour factors proportional to the cubic Casimir operator 
d abc d abc . Individually these diagrams are non-zero but this higher qroup invariant cancels in 
the sum. 

To obtain the 4-loop contribution to the anomalous dimension 7 s we need to calculate 
the MS renormalization factor Z^ g of Eq. ( |ll|) (and of Eq. ([Dj) where the lower order Zf^a 
appeared) which contains apart from a leading constant 1 only poles in the regularization 
parameter e, 

Z^ s (a s , e) = 1 + Z^ l \a s )/e + Z^ 2 \a s )/e 2 + ■■■ (15) 

The anomalous dimension 7 s is then expressed through the coefficient in front of the first pole 
in ^m S 

T s = -a s (j^Z^\a s )\ (16) 



To obtain the renormalization factor Z|^ s we calculated the overall ultraviolet divergence of 
the Green function 



G 



tl4]i> 



dxdye^ +i Py(0\T{^(x)J 5 JyM0)}\0) 



(17) 



This Green function is renormalized multiplicatively and the standard MS renormalization 
factor (containing only poles in e) is equal to Zf AS /Z 2 , where Z 2 is the renormalization factor 
of the inverted quark propagator. 

To obtain the required anomalous dimension 7 s in the af order we calculated the overall 
divergences of both the quark propagator and the Green function G^r j5 i , in the 4-loop order. 

This can be conveniently done using the technique that is described in Ref. [^8| . This general 
technique is based on the direct calculation of 4-loop massive vacuum (bubble) integrals and 
provides a procedure that is well suited for the automatic evaluation of huge numbers of 
Feynman diagrams. For the present calculation we needed to evaluate of the order of 10000 
4-loop diagrams. These diagrams were generated with the program QGRAF |p5f|. 
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27 
560 s 

~8T' 



3 -( 3 ) + C A C F nfT F (- m 



58 

y 
+ 



368 
-3-C3) 
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(18) 



The results of Eqs,( |i4| , |l8| ) are obtained in an arbitrary covariant gauge for the gluon field. 
This means that we keep the gauge parameter £ that appears in the gluon propagator i 
[—g^ + (1 — CjQ^tf /(q 2 + ^ e )]/('7 2 + ^ e ) as a f ree parameter in the calculations. The explicit 
cancellation of the gauge dependence in the coefficient function and the anomalous dimension 
gives an important check of the results. The results for individual diagrams that contribute to 
C s and 7 s also contain (apart from the constant £3) the constants £4, £5- The cancellation of 
these constants at various stages in the calculation provides additional checks of the results. 



We should also note that the various higher order colour factors [see Ref. |28]] that appear 
in the separate results for Z 2 and (Z^g/^) all canceled in Eq. (JT8J) . 

At this point it is interesting to compare the obtained result for 7 s with the analogously 
defined non-singlet anomalous dimension 7 ns since this can be constructed at 4-loops from 



known 3-loop results. Since the non-singlet anomalous dimension 7 ns vanishes, we have [see, 
e.g. Eq. (|l^)1 7 ns = —f3{a s )-^-log{Z^') where Z£ s is the finite renormalization constant for 
the non-singlet axial current which is known in the order a^. Since the difference between 
the singlet and non-singlet sector is in terms of at least one power of nj we must have that 
all the n^-independent terms in 7 ns and 7 s coincide, and indeed they do. This gives another 
strong check to the calculations. 

Substitution of the obtained anomalous dimension and coefficient function in Eq. (|9j) gives 
the following result for the Ellis-Jaffe sum rule. 
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(19) 



where a s = a s (Q 2 ) = 4-7ra s (Q 2 ), /3o = 11 — 2/3nj is the 1-loop coefficient of the beta function 
and ao is the invariant matrix element of the singlet axial current defined in Eq. ([!(]) . 
In particular, for nj = 3 we find 



dxg p l w {x ) Q' 2 
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3.5833 



ru 



IT 
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ru 



IT 



{± T2 l9Al + h a8) 
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IT 



0.54959 



a. 



7T 
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a (20) 



It is interesting to compare our result Eq. ( [19| ) with a recent estimate J3(|] of the singlet 
coefficient d| for n^ = 3. The estimate gives the value -2 which is slightly less than half the 
value obtained in the present article. 

In table 1 we have listed the numerical values of the second and third-order coefficients 
for the Ellis-Jaffe sum rule for n/ = 3, 4, 5, 6. One can observe the sign-constant charac- 
ter of perturbative QCD series both for non-singlet and singlet contributions. The series 
tends to preserve its sign-constant character even when perturbative coefficients of the singlet 
contribution change their signs around the value nj = 4. 



n f 


non- 


-singlet 


singlet 






KAO 2 


(a s /7r) 3 


(«» 2 




(a s /n)* 


3 


-3.58333 


-20.21527 


-0.54959 




-4.44725 


4 


-3.25000 


-13.85026 


1.08153 




4.87423 


5 


-2.91667 


-7.84019 


2.97845 




13.07103 


6 


-2.58333 


-2.18506 


5.27932 




20.73034 



Table 1. Second and third-order coefficients for the Ellis- Jaffe sum rule. 



Another deep inelastic sum rule, the Bjorken sum rule for neutrino-nucleon scattering, 
which is also known in the a^-order [37], also exibits the sign-constant behaviour of the 
perturbative QCD series. One can see that the obtained perturbative coefficients of the Ellis- 
Jaffe sum rule grow rather moderately. If we assume that the error of the truncated asymptotic 
series is determined by the last calculated term, then the obtained af approximation for this 
sum rule provides a good theoretical framework for extraction of the fundamental constant 
a,Q from experiment. 
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